Electronic Correlations in Oligo-acene and -thiophene Organic Molecular Crystals 
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From first principles calculations we determine the Coulomb interaction between two holes on 
oligo-acene and -thiophene molecules in a crystal, as a function of the oligomer length. The elec- 
tronic polarization of the molecules that surround the charged oligomer, reduces the bare Coulomb 
repulsion between the holes by approximately a factor of two. The effects of relaxing the molecular 
geometry in the presence of holes is found to be significantly smaller. In all cases the effective hole- 
hole repulsion is much larger than the valence bandwidth, which implies that at high doping levels 
the properties of these organic semiconductors are determined by electron-electron correlations. 

PACS numbers: 71.20.Rv,71.27.+a,71.38.-k 



A unique aspect of organic molecular crystals is that 
the electronic properties of these solids bear the marks 
of both molecular and condensed matter physics. In 
a crystal the organic molecules preserve their identity, 
since molecular crystals are held together by van der 
Waals interactions without the formation of intermolec- 
ular covalent bonds. This results in a small overlap be- 
tween the wave functions of electrons located on adjacent 
molecules, which leads to rather narrow electronic bands. 
Because of the small electronic bandwidths in these ma- 
terials, the energy associated with interactions such as 
the electron-phonon or electron-electron interaction, may 
dominate over the kinetic energy of the charge carriers. 

At a low density of charge carriers the interaction with 
phonons is expected to be dominant, which possibly leads 
to the formation of polarons At increased carrier 
density -approaching one carrier per molecule- the role 
of electron-electron interactions becomes more and more 
important. In this regime Coulomb interactions may sta- 
bilize correlated magnetic ground states. 

Probably the most studied example of organic conduc- 
tors in the high-density regime is that of Ceo ■ In this ma- 
terial it has been shown experimentally that the Coulomb 
repulsion between two carriers on the same molecule is 
substantially larger than the electronic bandwidth 0. 
At high carrier concentrations, the large Coulomb inter- 
action results in a strongly correlated ground state. Both 
K 4 C 60 and ammoniated K 3 C 6 o are Mott-Hubbard insu- 
lators, for instance 0,0]. 

It can be expected that electron-electron interactions 
are not only important in Cgo crystals, but also in crys- 
tals of many other molecules. In this paper we con- 
sider oligo-acenes and oligo-thiophenes, two classes of 
molecules often used in organic electronic applications, 
and we study their electronic structure in the high den- 
sity regime using first principles calculations ||. 

For all the molecules of experimental interest we find 
that the Coulomb interaction between two carriers on 



FIG. 1: Oligo-acenes are ladder-type hydrocarbons and 
oligo-thiophenes consist of a chain of thiophene rings. The 
molecules (a) tetracene and (b) quater-thiophene are shown, 
which have n = 4 rings. They are planar and have (a) Din 
and (b) C-zh symmetry. The black, grey and white balls rep- 
resent sulphur, carbon and hydrogen atoms. 



the same molecule is much larger than the electronic 
bandwidth, also when the electronic screening of this in- 
teraction in the crystal is taken into account. Our re- 
sults indicate that, analogous to C 60 , crystals of oligo- 
acenes and -thiophenes should be viewed as correlated 
electron systems. In addition we find that, compared to 
the Coulomb interaction, the energy associated with re- 
laxing the molecular geometry upon charging is smaller. 
This implies that the charge-charge Coulomb repulsion 
is the dominant factor at high carrier density. Correla- 
tion effects in these systems should be even more pro- 
nounced than in Cgo ■ Many of the electronic molecular 
levels in Ceo are degenerate due to the high symmetry 
of this molecule and such degeneracies effectively reduce 
the Coulomb interaction In contrast, the symmetry of 
oligo-acene and -thiophene molecules (see Fig.^) is much 
lower and all their molecular levels are non-degenerate. 

Bare molecular Coulomb interactions. We first present 
the results obtained for the Coulomb interaction strength 
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[/bare between holes on single, isolated oligo-acene or - 
thiophene molecules of different sizes. They are indicated 
by n, the number of (benzene or thiophene) rings com- 
prising the oligomer, cf. Fig. 2] To establish the value 
of t/bare we nrs t consider two singly charged n + ions far 
apart. The bare Coulomb interaction is then determined 
by the process in which a single electron is moved from 
one molecule to the other, i.e. n + + n + — ► n + n 2+ . 
Consequently the molecular Coulomb repulsion for two 
holes can be expressed in terms of the total energies 
of the neutral, doubly and singly charged molecules as 
[/bare = E(n) + E{n 2 +) - 2E(n+). 

Only ground state energies of individual molecules and 
ions enter here, which can be obtained with high ac- 
curacy from first-principles electronic structure calcula- 
tions. The results presented here are obtained by density 
functional theory (DFT) in the local density approxima- 
tion; using a generalized gradient functional gave only 
small differences Q. Pseudo-potentials are used to rep- 
resent the atomic cores || and the valence electronic wave 
functions are expanded in a plane wave basis set. Molecu- 
lar calculations are performed by enclosing the molecule 
in a finite box 0, whereas the crystal calculations use 
periodic boundary conditions. This scheme has proved 
to yield accurate results for molecular crystals 01 , poly- 
thiophene and oligo-thiophene dimers [l^ . 

The calculated bare Coulomb interaction [/bare be- 
tween two holes on a n-acene or n-thiophene molecule, 
n ranging from 2 to 5 and from 4 to 16, respectively, is 
listed in Table I, where the molecules are ordered accord- 
ing to increasing length, [/bare ranges from 5.3 to 1.5 eV; 
it decreases as the length of the molecule increases. On 
conjugated molecules the holes are delocalized over the 
molecule, so their average distance is larger on a larger 
molecule, which of course lowers their Coulomb repulsion. 

Molecular relaxations. So far we have not incorporated 
the fact that a charged molecule relaxes its geometry, 
which lowers its energy. This relaxation depends on the 
charge state of the molecule and effectively reduces the 
bare Coulomb interaction. The relaxation energy is the 
energy gained when optimizing the geometry of a charged 
molecule, starting from the geometry in its neutral state. 

As for the bare Coulomb interactions, the relaxation 
energies also decrease with increasing molecular length n. 
Upon charging an oligo-thiophene molecule, it changes 
from an aromatic to a quinoid geometry [lOj. A good 
measure for characterizing the change is to monitor 
the difference Ar between the C1-C2 and C2-C2 bond 
lengths, cf. Fig. ^ where Ar < for the aromatic struc- 
ture, and Ar > for the quinoid structure. The decrease 
of the relaxation energy -Ereiax with increasing length n 
can be understood qualitatively assuming that the ex- 
cess charge is delocalized over the molecule. The induced 
change of Ar should then roughly be proportional to the 
average excess charge density, i.e. Ar oc Z/n, where 
Z = 1 [Z = 2) for the singly (doubly) charged molecule. 





n-acenes 


n-thiophenes 


n 


2 3 4 5 


4 6 8 12 16 


f^relax 
i-T AMI 
^ relax 

f/eff 


5.26 4.51 4.01 3.63 
5.10 4.42 3.95 3.59 
5.00 4.32 3.86 3.52 
2.47 1.99 1.75 1.59 
1.22 1.23 1.15 1.10 


3.54 2.82 2.46 1.79 1.48 
3.34 2.71 2.37 1.69 1.38 
3.05 2.25 1.85 1.79 1.48 
1.47 1.21 1.09 0.70 0.56 
0.95 0.80 0.69 0.52 0.44 


u eB 

W 


2.67 1.96 1.65 1.39 
0.35 0.34 0.39 0.61 


1.43 1.10 0.98 0.64 0.50 
0.36 0.34 0.31 0.30 0.30 



TABLE I: Coulomb interaction U between two holes on n- 
acene and n- thiophene molecules, [/bare denotes the bare 
interaction; Kciax and [/^Ikx denote the relaxed interaction 
calculated with DFT and AMI, respectively (see text, H3)j 
L/fff, and C/ e ff denote effective interaction, calculated with the 
continuum model and the discrete dipole model, respectively. 
E^ ol denote the polarization energies calculated from the dis- 
crete dipole model. The bottom row lists the valence band- 
widths of the molecular crystals. All energies are in eV. 

Fig. |21 in which Ar obtained from the optimized geom- 
etry is plotted against Z/n, demonstrates that the rela- 
tionship is indeed approximately linear. In a harmonic 
approximation the relaxation energy per thiophene ring 
is proportional to (Ar) 2 . For an oligomer consisting of 
n rings, we then find E^ c ^ K (n) oc n(Ar) 2 oc Z 2 /n. The 
calculated relaxation energies indeed follow this scaling. 
Similar geometry changes take place upon charging the 
oligo-acene molecules, although the pattern is less uni- 
form. Their relaxation energies are much smaller, but 
also follow the EfJ/ {n) oc Z 2 /n behavior. 

Yet the absolute value of the relaxation energies is 
small compared to the Coulomb interaction. The "re- 
laxed" Coulomb interaction is defined as [/ re iax = [/bare — 
(iVelax — 2-^rciax) ■ The inclusion of molecular relaxation 
upon charging changes the Coulomb interaction by less 
than 10%; compare [/bare and [/relax in Table I. 

Meanwhile, it has been argued that DFT calculations 
underestimate relaxation energies |l3|. Therefore, we 
have also calculated the latter using the AMI semi- 
empirical Hartree-Fock method 0], which is known to 
yield a substantially larger relaxation and is likely to give 
an upperbound to the relaxation energies. The inclusion 
of AMI instead of DFT relaxation energies leads to !/M 
in Table I. For oligo-thiophenes the AMI relaxation ef- 
fects are indeed much larger 0, , whereas for oligo- 
acenes these effects are marginal. Note however that, ir- 
respective of the technique used to calculate relaxations, 
the Coulomb interaction is dominant. 

We therefore conclude that in the high density regime 
the nature of the ground state and low energy quasi- 
particle excitations is foremost determined by electron- 
electron interaction effects. The electron-phonon cou- 
pling is then a secondary effect which may result in a 
renormalization of quasi-particle properties. 
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FIG. 2: The bond length change upon charging of n-thiophene 
molecules as function of the average excess charge density 
Z/n. The open and closed circles represent data from singly 
and doubly charged molecules; the line is a linear fit. 



Screened Coulomb interactions. Having discussed the 
response of a single molecule to the creation of holes, we 
will now consider the response of the molecular crystal. 
In the presence of excess charges the electronic system of 
the crystal reacts to effectively screen the Coulomb inter- 
action between the charges. The electronic screening due 
to electrons on the same molecule on which the charge(s) 
reside is already implicitly accounted for in the calcula- 
tion of J/barc- We determine the screening energy due to 
all molecules that surround a charged molecule by cal- 
culating the polarization of these molecules in a discrete 
dipole model, see Refs. HillrR If we denote the total po- 
larization energy that results from singly/doubly charged 
molecules by E , , the effective screened Coulomb in- 

J pol 3 

teraction is given by U cS = U rciax - (E^+ - 2£+ ol ). 

I /2 I 

In order to evaluate the polarization energy E p ^ , 
one needs to sum over the electric dipole moments on 
all molecules of the crystal lattice, induced by the excess 
charge 0, 0] . The strength of an induced dipole is pro- 
portional to the local electric field, which is the sum of 
the field due to the charge and of that due to all the other 
induced dipoles. The proportionality factor is the molec- 
ular polarizability tensor, whose values for oligo-acenes 
and -thiophenes are taken from Refs. and llfll The 
lattice sum is determined by the crystal structures, which 
are taken from Refs. r2(ir2lU'2'2t r2.'j[ Large molecules are 
not well represented by point charges or point polariz- 
ability tensors. We therefore use the sub-molecular ap- 
proach, where the charge or polarizability is distributed 
over the molecule by assigning a (sub-molecular) charge 
or polarizability to each benzene or thiophene ring [Hll8|. 

One complicating factor for 4- and 6-thiophene is that 
at room temperature two crystal structures exist, called 
the LT and HT phases In practice the calculated 

E* oi turns out to be the same within 0.01 eV for these 
two phases, so we will not list them separately. Experi- 
mental crystal structures for 12- and 16-thiophene are not 
available and, for purpose of comparison only, we artifi- 
cially construct their crystal structures by starting from 



the HT phase of 6-thiophene and enlarging the crystal 
axis that points along the long molecular axis, such as to 
accommodate the longer 12- or 16-thiophene molecule. 

The results for the polarization energy E pol are shown 
in Table I. The values are between 1.2 and 0.4 eV, 
where increasing the length n of the molecule decreases 
the polarization energy. This is basically a volume ef- 
fect; for longer molecules the charge is distributed over a 
larger volume, which "dilutes" its electric field. Whereas 
the polarizability per unit volume of the surrounding 
molecules increases somewhat with the size of the conju- 
gated molecule, this dependence is weak and it does not 
compensate for the weaker average electric field. 

The polarization energy is slightly different if one uses 
the ionic charge distribution given by the AMI calcula- 
tions. This is more localized than the DFT charge distri- 
bution, leading to a higher E^ ol . In practice we observe 
this only in the larger thiophenes, where E^ ol becomes 
0.87, 0.79 and 0.72 eV for 6-, 8- and 12-thiophene. 

At this point we are able to determine the effective 
Coulomb interaction U e ff and the results are shown in 
Table I. Using AMI charges, U c s changes by less than 
0.06 eV even for the larger thiophenes, so we have not 
listed these AMI results. One observes that screening 
due to polarization of the surrounding molecules has an 
effect on the Coulomb interaction which is almost an or- 
der of magnitude larger than that of molecular relaxation. 
In fact, the screening is such, that the effective Coulomb 
interaction U c s is less than half [/bare- 

Since polarization has such a large effect, one would 
like to establish the sensitivity of the results upon the 
model used to calculate the polarization energy. There- 
fore we have also calculated the polarization energies 

j In j 

J5 j using a conceptually much simpler, but cruder 
continuum model. In such a model the charged molecule 
is placed inside a cavity surrounded by a homogeneous 
dielectric medium, where the cavity is shaped accord- 
ing to the molecule. The polarization energy is then 
given by the electrostatic interaction between the molecu- 
lar charge distribution and the surrounding medium |l2| . 
Table I also contains the results for the screened interac- 
tion U^g based upon polarization energies calculated us- 
ing a homogeneous dielectric medium. We have used a di- 
electric constant e = 3 for the oligo-thiophenes, as well as 
for naphthalene (2-acene) and e = 3.3 for the other oligo- 
acenes 18] . One observes that the screened interaction 
calculated with the crude continuum model is roughly 
within 10% of the microscopic discrete dipole model. 
Moreover, U ha , TC - U^ s = E^ - 2E+ ol scales as (e - l)/ e 
so it is not very sensitive to moderate changes of e [l2| • 
We conclude that the values of the screened Coulomb in- 
teraction obtained by the discrete dipole model are not 
very sensitive to the details of that model. 

Coulomb interaction versus bandwidth. To determine 
the correlation strengths, the effective Coulomb inter- 
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20 40 
length (A) 

FIG. 3: The effective Coulomb interaction t/ e ff and valence 
bandwidth W (Wlt for the LT phase) as function of the 
oligomer length (defined as the distance between the two ex- 
tremal carbon atoms). The data for the oligo-acenes and - 
thiophenes are plotted left and right of the solid vertical line. 



action U c s for holes should be compared to the band- 
width W of the highest valence band. The latter can 
be obtained from DFT calculations on oligo-acene and 
-thiophene crystals. Crystals of the HT and LT phases 
of 4- and 6-thiophene have been produced, whereas 8- 
thiophene is only reported in the LT phase Again, 
for purpose of comparison only, the crystal structures 
for the HT phases of 8-, 12- and 16-thiophcnc arc artifi- 
cially constructed in the way discussed above. The cal- 
culated valence band widths W are shown in Table I. Val- 
ues in the same range have been reported on the basis of 
semi-empirical calculations [3 US • Generally, the band- 
widths depend only weakly on the size of the oligomer. 
For instance, it decreases from 0.36 eV to 0.30 eV going 
from 4- to 16-thiophcne in the HT phase. Such a weak 
dependence can be explained by the observation that the 
highest valence band is mostly derived from the highest 
occupied molecular orbital and the latter has the same 
character for all oligo-thiophcncs. The bandwidth in the 
LT phase (not shown in Table I) decreases from 0.22 eV, 
0.21 to 0.16 eV for 4-, 6- and 8-thiophcne, respectively. 
These values are smaller than those of the HT phase be- 
cause of the different packing of the molecules in the two 
phases. The bandwidths of the oligo-acenes are similar 
to those in the oligo-thiophenes; only that of pentacene 
is substantially larger. 

Fig. |21 summarizes the main result of our study. The 
figure shows the calculated effective Coulomb interaction 
Z7 e ff and the bandwidth W for the oligo- acenes and - 



thiophenes, as a function of their molecular lengths. For 
longer molecules U e s becomes smaller, but does not drop 
below W for any size of molecule. For the molecules that 
are of main experimental interest at present, i.e. 4-, 6- 
thiophene and the acenes, the ratio U e s/W is between 
2.3 and 7.6. So we conclude that in these systems the in- 
tramolecular Coulomb interaction is very large compared 
to the electronic bandwidth. We therefore expect that 
when oligo-acene and thiophene molecular crystals are 
highly doped, correlation effects will be very prominent, 
giving rise to magnetically ordered states. 
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